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Abstract 

We explore constraints on the spectral index n of density fluctua- 
tions and the neutrino energy density fraction 0,hdm, employing data 
from a variety of large scale observations. The best fits occur for n ~ 1 
and 0,HDM ~ 0.15 — 0.30, over a range of Hubble constants 40 — 60 
km s~^ Mpc^^. We present a new class of inflationary models based 
on realistic supersymmetric grand unified theories which do not have 
the usual 'fine tuning' problems. The amplitude of primordial density 
fluctuations, in particular, is found to be proportional to (Mx/Mp)^, 
where Mx{Mp) denote the GUT (Planck) scale, which is reminiscent 
of cosmic strings! The spectral index n = 0.98, in excellent agreement 
with the observations provided the dark matter is a mixture of 'cold' 
and 'hot' components. 

PACS Nos: 98.80Cq, 98.65Dx, 12.10Dm, ll.SOPb 
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Recent studies of large scale structure formation [|I|, when confronted with 
a variety of data from the Cosmic Background Explorer @ (COBE) and other 
large scale galaxy surveys, provide support for an inflationary scenario in 
which the spectral index of density fluctuations n is close to unity and the 
dark matter is a mixture of cold and hot components. Non-supersymmetric 
grand unified theories which give rise to precisely this scenario were con- 
structed more than a decade ago [Q. However, several fundamental chal- 
lenges, including unification of GUTS with gravity and the gauge hierarchy 
problem, strongly hint that the supersymmetric grand unified (SUSY GUTS) 
framework may be a more promising way to proceed. 

Supersymmetric GUTS based on G = SU{?>)c x SU{?>) l x SU {?>)r have the 
desirable feature that they permit unification of the standard model gauge 
couplings to occur at scales on the order of IQ^^GeV , which is indicated by the 
recent LEP data. Moreover, as has recently been shown j^, by introducing 
either some symmetries, or alternatively R-symmetry, the gauge hierarchy 
('fine tuning') problem can be resolved in SUSY GUTS based on G. It should 
be stressed that this is accomplished without complicating the 'higgs sector' 
of the theory. The 'minimal' choice will do, which is well nigh impossible 
in GUTS such as SU{b) or 50(10). The additional symmetries also ensure 
that the proton is essentially stable. 

Encouraged by these developments, we investigate here if the inflationary 
scenario can be realized within the framework of G. We would call the 
attempt 'successful' if the following conditions are met. First, the scalar 
(higgs) sector of the theory, including the infiaton part, is determined by 
particle physics considerations. Second, no 'fine tuning' of parameters is 
needed. [This includes the gauge hierarchy problem.] Finally, the scenario 
must be such that the Planck scale corrections from supergravity and/or 
superstrings can be safely ignored. Although non-trivial, it turns out that 
all of these constraints can be satisfied within the framework of G. The 
inflationary scenario we are led to has previously been considered in general 
terms by Linde [^] (see also ref. [^), and been dubbed 'hybrid' inflation. 
Our realization of 'hybrid inflation' within a supersymmetric framework is 
unique in a number of ways and can be implemented in a variety of models, 
especially those based on G. One particularly important result has to do 
with the amplitude of primordial density fluctuations, which turns out to be 
proportional to (M^/Mp)^, where Mx denotes a superheavy (GUT scale) 
and Mp ~ 1.2 x 10^^ GeV is the Planck mass. The spectral index of the 
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density fluctuations is very close to unity as required by the observations. 

We begin with a definition of the cosmological parameters in the simplest 
inflationary models. 

1. The density p is the critical density pc = ?>Hll {9>tiG), where Hq = lOOh 
km s~^ Mpc~^ is the present value of the Hubble constant and G is 
Newton's gravitational constant. 

2. We take h = 0.5 ± 0.1, which corresponds to the range of h allowed by 
the combination of observations and constraints on the age of the uni- 
verse in critical density models with a vanishing cosmological constant. 

3. We take the baryon fraction to be the midpoint of the range allowed 
by big bang nucleosynthesis 0, flbaryonh"^ = 0.0125 ± 0.0025, where 

^baryon — Pbaryon/ Pc- 

4. The remaining mass density is in the dark matter which is composed 
of some mixture of cold (CDM) and hot (HDM) components. The 
latter is assumed to be "lightly" massive (few eV) relic neutrinos. The 
relative concentrations of the two components is left to be determined 
by the data. We include the possibility that the dark matter may be 
all CDM. 

5. The power spectrum, which is the Fourier transform of the two-point 
density autocorrelation function, has the primordial form P{k) oc /c", 
where k is the amplitude of the Fourier wavevector, and n denotes the 
spectral index. The value n = 1 corresponds to the Harrison- Zeldovich 
spectrum. 

The above set of parameters deflne a family of inflationary models. Their 
ultimate test comes from comparison to data. Here we will extend the treat- 
ment previously described in Ref. [|I| to include larger values of the density 
fluctuation spectral index n and different values of the Hubble constant. We 
will give a brief summary of the procedure used there and present the new 
results in Figure 1. The data we use is summarized next. 

We begin with the data from COBE. The amplitude is most conveniently 
characterized in a relatively n independent way by the expected hexade- 
cupole moment of the temperature anisotropy 6T4 = 12.8 ± 2.3 pK, where 
the error includes effects such as radiometer noise and cosmic variance [Q]. 
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(This treatment is consistent with the recent analysis of the COBE two year 
data 0.) This data set probes the spectrum on the largest scales (~ 10'^ 
Mpc). For smaller scales, we use the data obtained from large galaxy sur- 
veys, particularly that of the Infrared Astronomy Satellite (IRAS) catalogue 



of galaxies. In ref. this catalogue was used to obtain a direct estimate 



of the power spectrum in the range 30 /i/Mpc < 27r/fc < 300 /i/Mpc. The 
resulting power spectrum is in good agreement with estimates from other 
surveys based on optically selected galaxies. The IRAS catalogue of galaxies 
has also been used in a different way to determine the large scale bulk velocity 
field. We use the values given by the potential flow algorithm ("POTENT"' 



collaboration (see ref. ||11|| ). 



There are important constraints on the power spectrum which come from 
structures which have evolved so that their description is inherently non- 
linear. If we make some basic assumptions about the formation of these 
structures, we can use this data to constrain the amplitudes of the fluctua- 
tions in the linear theory. We will use two such pieces of information. First, 
in pre-COBE days the normalization of power spectra was usually done by 
specifying the mass fluctuation on 8 Mpc scales, and much effort has 
been made to determine this normalization from the data. In ref [12|, an 
attempt was made to estimate this quantity from the abundance of clusters. 
Estimates of the mass fluctuation on 8 Mpc from different galaxy sur- 



veys yield results that are consistent [|^. We adopt the synthesized value 



SM/M{8 Mpc) < 0.8 as a constraint on the amplitude of linear pertur- 
bations. 

Second, the number of early quasars at high redshift requires a minimum 
amplitude of mass fluctuation. This issue has recently been extensively stud- 
ied Jl^. In Ref. two of the present authors estimated the error on this 
amplitude introduced by theoretical modeling uncertainties and found that 
the lower limit on the amplitude of mass fluctuations required to make the 
early quasar population is 6M/M{0.6 Mpc) > 1.1 ±0.2. We will use this 
constraint here as well ||15|| . 

The data are then analyzed in the following way. For a given value of 
n and a hot dark matter fraction Qhdm, we calculate up to a normaliza- 
tion factor, the hexadecupole moment of relic temperature fluctuations, the 
power spectrum for the IRAS power spectrum wavenumbers, the large scale 
velocities, the linear mass fluctuation on 8 Mpc, and the amplitude of 
the quasar sized mass fluctuations. In addition, while most of our list of 
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constraints depend on the mass fluctuations, the variations in galactic num- 
ber can be somewhat different than mass fluctuations. To compare our mass 
fluctuation power spectrum with galactic number we have to also introduce 
a "density bias factor" bj to allow for the possibility the IRAS galaxies may 
not trace the mass. The fact that the depth coordinate in the galaxy survey 
is doppler redshift space rather than physical length introduces some effects 



which must be corrected for (see ref. [|T6l). For each value of VLhdm and n 



we first do a least squares fit of the normalization and the bias factor 6/. 
After fitting these parameters, we proceed to calculate the value of the 
(goodness of fit) statistic from the data for each model. For the bulk 
velocities, we have combined the uncertainty due to "cosmic variance" in 
quadrature with the quoted analysis errors, because the velocities correspond 
to a very limited sampling of a Gaussian velocity field. After finding the 
best fitting model, we record how the statistic changes as we vary the 
parameters n and fin dm and draw confidence level contours in the ti — VLhdm 
plane [l^ . We repeat this entire procedure for different values of the Hubble 
constant. 

Figure one dramatically shows how the data favor values of n very close 
to unity with the actual limits depending somewhat on the mix of dark mat- 
ter and the Hubble constant. Overall we can say that with 99% confidence, 
0.80 < n[i^o/(50km s~^ Mpc"^)]^''^ < 1.15, independent of the dark matter 
composition. The value of n = 1 works well over the range of Hubble con- 
stants 40 — 60 km s~^ Mpc~^ It appears that n is constrained to be quite 
close to unity by the data. 

We now discuss how such density fluctuations can be realized in realistic 
super symmetric GUTS. We are particularly interested in identifying models 
in which there are no 'fine tuning' (including gauge hierarchy) problems. To 
set things up, consider the following globally supersymmetric renormalizable 
superpotential W: 

W = K S ^ (j)- ^i^S (1) 

where 0(0) denote a conjugate pair of superfields transforming as non-trivial 
representations of some gauge group, while is a gauge singlet superfield. 



This superpotential is 'natural' in the strong sense |18[. It is of the most 
general form consistent with R-symmetry under which 5 — > e*'^S', — *• 
e^'^W, while the product 00 is invariant. Note that cubic terms in and can 
be forbidden by assuming, for example, the transformations —>■ e*'''0, — *• 
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e^"'(j). In realistic models such terms may be allowed without altering the 
main conclusions. 

We point out that, at least in the global SUSY case, the R-symmetry is 
the unique choice for implementing the 'false' vacuum inflationary scenario 
in a natural way. It is the only symmetry which can eliminate all of the 
undesirable self-couplings of the (inflaton) S*, while allowing the linear term 
in the superpotential. With supersymmetry unbroken, the potential takes 
the form (we represent the scalar components with the same symbols as the 
superfields if there is no danger of confusion!): 

+ I - /i^ P +D - terms (2) 

The D-terms vanish along the (D-flat) direction | | = | 0* |. Consequently, 
the only supersymmetric minimum of the potential is at 

< ^ > = 

(3) 

Mx = {\<p\) = {\^\) = /i/v^(^>0,«:>0) 

We will say more about the scale Mx shortly. 

Consider now an early universe scenario with chaotic initial conditions. 
For I S \>\ Sc 1= yu/v^) the effective potential V is minimized by 
< (f) >=< 4> >= 0. That is, for 15*1 > Sc, the energy density is dominated 
by the 'false' vacuum energy density /i^, which can therefore lead to an ex- 
ponentially expanding (inflationary) universe. The potential in (2) does not 
contain a term which can drive S to its minimum value. This, however, is no 
longer the case when the quantum corrections are taken into consideration. 

With l^l > Sc, both and vanish and there is a nonzero Fs — term{= 
/i^) which breaks supersymmetry, such that the one loop corrections to the 
effective potential are non-vanishing, and given by []TU[ 

where the summation is over all helicity states, (—1)^ indicates that the 
bosons and fermions make opposite sign contributions, and A denotes a renor- 
malization mass. The quantum corrections will help drive S to its minimum. 

Note that for S > Sc there is no mass splitting inside the gauge supermul- 
tiplets or the S-superfield (actually the masses of S-scalar and its fermionic 



V{S, 



s 



6 



superpartner both vanish). The non- vanishing contribution is from the mass 
sphtting within the (f) , (f) superfields. The complex scalars in 0, are spht 
by the nonzero F^-term into two pairs of real scalar and pseudoscalar com- 
ponents with mass squared k'^S'^ ± k/x^, whereas the fermionic partners have 
mass kS. The one loop corrected effective potential (along the inflationary 
trajectory S > Sc, (j) = 4> = 0) is given by 

{kS' - f^'flnil - ^) + i'^S' + f^rinil + ^)] (5) 
If S is sufficiently greater than Sc, Veff{S) reduces to the simpler form 



A ^ S 3 . 

' + 32^5""^ + 2' 



(6) 



For 15"! > Sc, the inflationary phase is dominated by the false vacuum 
energy //^ as in the tree level case, but the additional contribution in (5) will 
now drive S to its minimum. The GUT phase transition takes place only 
after the S field drops to its critical value Sc{= Mx)- Below the S field is 
driven to zero by the positive mass term K^IS'pl^^l which is increasingly more 
effective due to the increase of the 0,0 vevs (induced by the decreasing S). 
All of the fields rapidly adjust to their vacuum values (3), thereby restoring 
supersymmetry. 

Note that the end of inflation does not necessarily coincide with the GUT 
phase transition which occurs when S approaches Sc- The end is signaled 
when the "slow roll" condition is violated for some S > Sc- We can charac- 
terize the "slow roll" condition as (see first paper in ref. [^) 



where 



e<<l,|?7|<<l (7) 



' ' V = ^TT (8) 



(the prime refers to derivatives with respect to S)- The inflationary phase 
may end before the GUT transition if the above conditions are violated at 
some S > Sc- For convenience, we can use the parametrization S = xSc, 
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where the parameter x characterizes the rolhng of S. (The GUT phase 
transition occurs for a; = 1.) The quantities e and r] are given by 



Note that rj becomes infinitely large for a; = 1, so that inflation ends as x 
approaches 1 (from above). 

So far, we have not introduced any supersymmetry violation in the system 
(the global minimum in (3) is supersymmetric). In conventional schemes 
(say N=l supergravity), this breaking is introduced through the soft SUSY 
violating terms in the tree level potential. The main influence of such terms 
on the inflationary scenario discussed above arises from the fact that the 
SUSY breaking induces a TeV scale (mass)^-term for the scalars, in particular 
for the S field. The term m^jS'p (m ~ TeV) provides an extra force driving 
S to the minimum. However, unless the coupling constant k is very small, 
the soft mass terms only provide a small correction to Veff{S) in (5), and 
so cannot significantly affect the above dynamics. This is not surprising 
since for \S\ > Sc, the non-supersymmetric {massY splitting inside the 0,0 
superfields is nfi^ which, as we shall see, is much larger then m^. In such 
a situation the inflationary scenario above is practically independent of the 
particular mechanism of supersymmetry breaking. 

Let us now compare the predicted quadrupole anisotropy, based on (5), 
with the value (^ 7 x 10^^) measured by COBE. From the scalar density 
fluctuations one has (see first paper in ref. [0]) 

\T )q ~ V 45 FMf I^Q (10) 

^ {MNqY^{Mx/MpY 

where the subscript xq indicates the value of S as the scale (which evolved 
to the present horizon size) crossed outside the de Sitter horizon during infla- 
tion, and A^q(~ 50 — 60) denotes the appropriate number of e-foldings. The 
formula in (10) is remarkable in that the fluctuation amplitude is propor- 
tional to (Mx/Mp)^, just as in the cosmic string scenario! The amplitude 
turns out to be in the right ball park, without having to 'fine tune' additional 
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parameters (such as dimensionless quartic couplings and/or the mass of the 
inflaton). Using (10), we can estimate the fundamental parameter Mx to be 
on the order of IQi^^-^GeV. We have ignored the contribution of the tensor 
fluctuations to the anisotropy, because they are suppressed by a factor of 
ft;/[87r(AiQ)^/^] relative to the scalar component in eq. (10). 
The spectral index n of the density fluctuations is given by 

n~l--^~0.98 (11) 

Nq 

which, as we have seen earlier, is in the central range of the values preferred 
by observations. 

An estimate of the coupling k is obtained from the relation 

A ^ (12) 
XQ 

With xq ~ 10 say, (which corresponds to S* ~ 10^^'^ GeV, and the Planck 
scale corrections can be safely ignored), the coupling k turns out to be on the 
order of 10"^. Note that for this value of n, the tensor generated anisotropics 
are less than 10~^ of the scalar anisotropy amphtude. 

Having outlined how super symmetric models can lead to a successful 
inflationary scenario without invoking small dimensionless couplings, we now 
proceed to discuss how this idea can be realized in realistic SUSY GUTS. 
We would like the model to be well motivated from the particle physics 
viewpoint. For instance, a non-trivial constraint would be the absence of 
the gauge hierarchy and proton decay problems. Although this constraint 
presumably can be met in GUTS such as SU{5) or 5*0(10), the higgs sector 
of the theory becomes quite messy. An even more formidable constraint 
stems from the fact that the phase transition involving the gauge non-singlet 
fields (0,0) occurs at the end of inflation. In SU{5) this would lead to the 
monopole problem. We therefore must resort to SUSY GUTS with rank 5 
or higher in order that at least the possibility exists for the monopoles to be 
inflated away. 

The simplest example of SUSY GUTS we are aware of in which the gauge 
hierarchy problem can be resolved with a minimal higgs system is based 
on G{= SU{3)c X SU{3)l x SU{3)r). Under the gauge group G, the left 
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handed lepton, quark and antiquark supcrficlds respectively transform as 
(1, 3, 3), (3, 3, 1) and (3, 1, 3). They are denoted as Aj, Qi and Qlii = 1, 2, 3): 

Ai — 



Qi 



d 1 (13) 



\9 ) 

Here H^^\ H^'^^ and L denote SU{2)l doublet superfields, e'^ is an SU{2)i 
singlet, and i/'^ are standard model singlets, while g{g^') denote additional 
down-type quark (antiquark) superfields. The SU{3)l {SU{3)r) symmetry 
acts along columns (rows). The symmetry breaking of G to SU{3)c x U{l)em 
requires at least two sets of higgs supermultiplets, A + A and A' + A', where 
A(A') transform as Aj above. The conjugate fields A(A') are needed to ensure 
that the SUSY breaking scale is well below the GUT scale and that the 
anomalies cancel. 

Imposition of a suitable R-symmetry, largely determined by the require- 
ment that a pair of electrowcak doublets (in A) remains light, fixes the renor- 
malizable part of the superpotential as follows (we will not give the couplings 
involving the chiral fermion families): 



W 



kSXX - iJ?S + 
6A'^ + cP 
ciA^A' + eAP 



(14) 



Here the bilinear and trilinear terms stand for the appropriate G invariant 
combinations. For instance, AA = Tr(A^A^), A'^ = e'^^^eABC>^'a>^f\'^, etc., 
where the Greek (Latin) indices refer to SU {3) l{SU {3) fc) ■ The R-charges of 
the superfields are as follows: ^"(1), A(— |), A(|), A'(|), A'(|), where W 
e^^W under R. The superlarge vevs will be oriented along the directions 



A Nl A* 




(15) 
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A' 1 = 1 X' 



/ 

I (16) 



V z/" 

The absence of a A'^ term in (14) guarantees that the electroweak doublet pair 
(H^^^ — H^'^^)m A remains massless at tree level, despite the GUT symmetry 
breaking. The well known '/i-term' for this pair is generated by the non- 
renormalizable couphng X^{\'\'Y/M^, (M = Mp/a/Stt ~ lO^^GeV^) which is 
the leading contribution allowed by the R-symmetry. With < A' > /M ~ 
10^^, fi is in the TeV range as desired. [Note that without the Wa,a' term, 
the superpotential in (14) would be invariant under [SU{3)l x SU{3)r]'^.] 

It is important to point out that in addition to resolving the gauge hierar- 
chy problem, the R-symmetry guarantees the absence of dimension four and 
five baryon number violating operators. That is, cubic superpotential cou- 
plings such as QiQjQk are absent, where the Qj's and Q"s carry R-charges 
= |. The proton is essentially stable according to this model. Moreover, 
the desired couphngs such as QiQ'jX, which provide masses for the quarks, 
are allowed. Similarly, there are no undesirable lepton number violating cou- 
plings. 

The superpotential (14) is the generalization of (1), on which the infla- 
tionary scenario was based, to a realisitic SUSY GUT model. Without going 
into the details, which will be published elsewhere, we briefly summarize how 
things work out. The role of the (j) — (j) — S sector in (1) is played by the 
superfields A — A — S*. [Note that the A^ term in (14) vanishes along the 
vacuum direction (15).] Prior to infiation, the vevs in the A' — A' sector, 
which arise from an interplay of suitable non-renormalizable couplings and 
supersymmetry breaking mass terms, break the {SU{3))^ GUT symmetry to 
SU (3)c X SU{2) L X SU{2y^ x U{1). The subsequent infiationary phase driven 
by the 5* field infiates away the monopoles produced during the first (A' — A') 
breaking. The A — A fields acquire non-zero vevs (~ 10^^'^ GeV) at the end 
of infiation. 

An infiationary scenario is incomplete without some discussion of i) how 
the baryon asymmetry is generated and ii) the dark matter candidate(s). As 
far as (i) is concerned, it has been shown in ref. [Q that the resolution of 
the gauge hierarchy problem within the framework of G necessarily leads to 
an additional pair of electroweak doublets. (Actually, one effectively ends 
up with a 5 + 5 of S'f/(5), of mass ~ TeV. Unification of the three gauge 
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couplings at scales ~ lO^^GeV is thereby preserved.) This strongly suggests 
generating the baryon asymmetry at the electroweak scale. The reheat tem- 
perature in this scheme is high, ~ 10^^ GeV, which means that the gravitino 
mass must be large in order that the primordial gravitino problem is avoided. 
Precisely how this works out would depend on how local supersymmetry is 
broken. 

As far as dark matter is concerned, the lightest supersymmctric particle 
(LSP) is completely stable and so presents itself as an ideal 'cold' dark matter 
candidate. Indeed, one can speculate here about the nature of the LSP. 
The point is that in the absence of the two lighter families, the well known 
parameter tanf3 (=< > / < H'^ >) of the MSSM is equal to mt/mh at 
the GUT scale Mx- (This follows from the superpotential coupling Q A.) 
Coupled with the radiative electroweak breaking scenario, this leads to a 
cold dark matter candidate essentially composed of the bino with mass in 
the range of 200-350 GeV. 

What about the neutrinos in this scheme? The lepton chiral supcrmul- 
tiplets = 1,2,3) transform according to the same G representation as 
the higgs supermultiplets A, A'. The R-symmetry allows the coupling AjAjA 
provided that the R-charge of A, is |. This is needed to provide mass say 
to the r lepton. But it also gives a Dirac mass m^tan/? to the (z/^ and z/^) 
neutrinos in Aj. Now, the R-symmetry prevents one from giving a large Ma- 
jorana mass to the fields, which could make the neutrinos unacceptably 
heavy. We are therefore led to introduce three G-singlet superfields Fj (one 
per chiral family) in order that the standard see-saw mechanism can be im- 
plemented. One of the neutrinos can now have a mass in the few electron 
volt as suggested by the data on large scale structure formation. 

In conclusion, it is gratifying to learn that a satisfactory inflationary 
scenario, consistent with large scale structure observations, can be realized 
in reahstic supersymmetric GUTS in which the sector that 'drives' inflation 
is determined by particle physics considerations, and no fine tuning of the 
parameters is needed. 
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1 Figure Captions 

Figure 1. The contours (68 %, 95%, and 99% confidence levels) for values 
of n and ^^hdm implied by the large scale structure data. In the left panel we 
have the constraints for a Hubble constant of 40 km s~^ Mpc~^. Similarly 
we have the constraints for Hubble constants of 50 and 60 km s~^ Mpc~^ in 
the center and right panels. Values of n ~ 1 are consistent with all values 
of the Hubble constant. We note that decreasing the Hubble constant favors 
larger values of n and smaller VLhdm- 
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